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a unique  motor  configuration,  different  from  existing  con- 
figurations, and  with  better  performance  characteristics. 

The  design  of  a generalized  electric  motor  is  formulated 
using  a doubly  excited  magnetic  field  system.  The  general- 
ized examples  of  transformer  and  synchronous  motor  designs 
lead  to  an  important  result:  there  exists  an  optimum 

machine  configuration  which  remains  unchanged  when  the  power 
of  the  motor  is  changed  over  a broad  range. 

The  design  of  a DC  PM  motor  is  optimized  for  minimum 
weight,  or  losses,  or  a combination  of  both  objective  func- 
tions. The  design  procedure  is  developed,  and  results  are 
compared  for  two  different  materials  (ferrite  vs.  rare 
earth  PM)  and  for  two  different  structures  (normal  vs. 
"inside-out") . The  comparison  shows  the  DC  PM  motor  with 
a rare  earth  magnet  in  an  "inside-out"  structure  with  a 
definite  "length/diameter”  ratio  to  be  a more  desirable 
solution  for  a robotic  application  because  the  weight  and 
inertia  are  reduced  noticeably  under  the  condition  of  reduced 
time  of  operation. 

The  time  of  operation  is  reduced  using  the  hierarchical 
design  procedure  and  using  the  constraint  of  maximum  bend- 
ing stress  of  the  gear  teeth.  The  time  of  operation  is  re- 
duced by  matching  load  and  actuator. 

The  results  of  the  research  help  to  improve  the  method 
of  electric  motor  design  as  well  as  the  performance  of  the 
mechanical  equipment  which  is  driven  by  the  motor. 


CHAPTER  I 
INTRODUCTION 

1 . 1 Problem  of  Optimum  Design  of  Electric  Motors 

Since  the  possibility  of  obtaining  mechanical  motion 
using  electromagnetic  means  was  discovered  by  Faraday  in 
1821  [11  , electric  motors  have  become  the  most  versatile  and 
efficient  devices  for  electromechanical  energy  conversion. 
Thus,  electric  motors  are  being  used  in  all  areas,  and 
design  techniques  have  been  developed  by  manufacturers  to 
improve  the  performance  characteristics  while  maintaining 
the  specifications  of  the  motors. 

The  performance  of  a motor  can  be  determined  by  its 
efficiency  [2-4],  acceleration  [5],  speed-torque  character- 
istics, and  the  temperature  of  the  winding  in  a variety  of 
work-regimes.  Most  of  these  parameters  are  related  to  the 
weight  and  losses  of  a motor.  Another  way  of  judging  the 
motor  performance  at  the  stage  of  motor  selection  among  many 
motors  is  to  compare  the  modified  parameters,  such  as  torque 
per  weight  (T/W)  [6]  or  power  per  weight  (P/W)  [7-11],  These 
parameters  can  better  characterize  the  energy  conversion 
capacities  of  the  motor  via  the  values  of  torque  or  power  to 
be  produced  by  the  unit  weight  of  the  motor.  The  specifica- 
tions of  the  motor  are  represented  by  output  power,  operating 
voltage,  and  rated  speed. 


selecting  principles 


The  word  “design”  means 
tion,  geometry,  dimensions,  and  materials  which  can  satisfy 
the  specifications  and  the  performance  requirements.  There 
is  an  infinite  number  of  possible  design  solutions.  The 

suitable  solution  among  the  multiplicity  of  admissible 
designs.  The  criteria  in  the  optimum  design  or  in  the  selec- 
tion of  a motor  from  the  list  of  available  choices  can  be 
oost,  weight,  or  losses  of  the  motor  combined  with  the  satis- 
faction of  the  requirements  of  the  work-regimes  (such  as  the 
temperature  or  the  time  of  operation) . Each  criterion  can 
be  used  as  a cost  function  or  an  objective  function  in  opti- 

The  approaches  to  the  design  procedure  changed  while 
the  status  of  motor  applications  changed  over  time  and  be- 
came more  diversified  in  its  purpose.  In  the  1920s  and  30s, 
one  motor  gave  mechanical  power  to  several  load  devices  using 
a belt  as  a power  transmission  device,  all  of  the  machines 

aged  the  individual  load  cycles  of  different  machines. 
Therefore,  at  this  time  each  load  was  not  to  be  specified. 

In  the  1930s  and  40s,  a single  motor  started  to  be  used  to 
handle  a single  machine;  so  users  began  to  individualize  the 
load  requirements.  There  were  still  many  different  mechani- 
cal operations;  so  it  was  necessary  to  average  the  load  with 
time.  The  ideas  of  “one-hour- load"  and  so  on  appeared  in 
the  practice  of  the  design,  simplifying  the  process  of 


dealing  with  a stochastic  load  on  the  shaft.  Beginning  in 
the  1950s,  several  motors  were  used  to  operate  a single 
machine  or  mechanical  system,  and  the  load  conditions  could 
be  determined  precisely  by  individualizing  each  motion  of 
this  system.  The  ultimate  example  of  this  system  is  an 
industrial  robot  where  sometimes  six  to  ten  motors  are  used 
in  the  robot,  and  each  motor  characterizes  each  motion  of 
the  system.  In  this  way,  the  applications  of  motors  have 
been  specialized  to  improve  the  motor  performance  by  matching 
each  load  and  motor  condition. 

The  motor  design  can  be  grouped  into  two  types  (as  shown 
in  Fig.  1-1),  which  correspond  to  the  status  of  knowledge  of 
the  individual  customer  and  his  requirements: 

TYPE  I:  Individual  customer  is  unknown  so  the  gener- 

al requirements  from  a group  of  similar 
customers  will  be  formulated. 

TYPE  Ila:  Individual  customer  is  known,  but  the  machine 

to  be  equipped  by  a motor  is  not  designed 

TYPE  lib:  Individual  customer  is  known,  and  the 

machine  to  be  equipped  by  a motor  is  already 
designed. 

In  TYPE  I design,  the  motor  will  be  independent  of  the 
customer's  requirements.  The  generalized  requirements, 
listed  below,  will  be  used  as  input  requirements: 

1.  The  average  load  torque  on  the  shaft  and  its  vari- 


2.  The  average  gear  ratio  in  the  machine  and  its 
variation. 

3.  The  average  inertia  of  the  machine  and  its  vari- 


4.  The  typical  cycles  of  operation. 

Under  all  these  conditions,  the  motor  should  be  designed  to 
be  cost  effective  while  maintaining  the  above  requirements. 


In  TYPE  Ila  design,  only  the  output  motion  of  the 
machine  to  be  designed  is  described,  and  the  load  character- 
istics are  given  as  follows: 

1.  (Load  torque,  speed)  = f (time) 

2.  Displacement 

3.  Velocity  trajectory 


In  this  case,  the  motor  torque  and  gear  ratio  are  selected 
to  provide  the  minimum  time  of  operation.  After  this,  the 
motor  should  be  designed  optimally  to  minimize  weight, 
losses,  or  cost  of  motor  under  the  obtained  motor  torque. 

In  TYPE  lib  design,  all  the  mechanical  information  about 
the  machine  can  be  obtained: 


1.  (Load  masses,  inertia)  - f (a,  1) 

2.  (Load  forces,  torque)  = f (a,  1) 

3.  Speed  and  dynamic  torque  dependencies  on  time  are 


This  motor  should  be  designed  to  minimize  weight,  losses, 
or  cost,  and  the  above  conditions  are  to  be  maintained. 


This  study  will  concentrate  on  the  motor  design  for 
actuators  in  robotic  manipulators  and  automated  manufactur- 
ing systems.  Since  the  robots  of  today  employ  almost  every 
kind  of  electric  motor,  the  technical  and  economic  character- 
ises robot  applications  produce  a clear  trend  toward 
higher  performance  drives. 

The  ideal  actuators  for  manipulators  must  provide  high 
productivity,  low  initial  cost,  accuracy,  uniformity  of 
motion,  low  losses,  and  high  reliability.  Actually,  it  is 
not  possible  to  achieve  the  best  values  of  all  of  these  prop- 
erties, so  each  superior  characteristic  requires  a compromise 
in  the  design  of  robotic  actuators.  The  economics  of  robot 
application  demand  that  a robot  have  a sufficient  productiv- 
ity to  pay  for  itself  in  a short  time.  The  actuator  systems 
in  present  manipulators  have  the  following  inherent  problems! 

1.  80-90%  of  work  is  for  carrying  its  own  weight. 

2.  40-50%  of  the  weight  is  in  actuators,  gearboxes 

and  supplemental  equipment. 

3.  There  is  a very  wide  range  of  dynamics. 

4.  The  main  factors  of  errors  are 

a.  Backlash  of  gears  . 

b.  Losses  of  energy  in  the  actuator  increases  the 
temperature  of  the  motor  structure  and  dimen- 
sional expansion  will  cause  performance  errors. 

c.  The  stiffness  of  the  actuator  system  will  have 
positional  errors  which  are  different  from  the 


The  theory  of  electric  motor  optimum  design  can  be 
tested  in  application  to  the  optimum  design  of  the  actuator 
for  a multilink  manipulator  to  solve  the  above  problems. 
Most  of  the  conditions  that  this  actuator  is  to  satisfy  are 
individual,  which  gives  us  an  opportunity  to  compare  the 
solution  with  several  conventional  alternatives  over  a 
variety  of  performance  indices  to  be  achieved. 

The  following  technical  requirements  can  be  set! 

!•  The  weight  of  the  motor  is  to  be  minimized  [11-121. 

2 . The  motor  is  to  provide  a variety  of  operational 
cycles  with  output  path,  speed,  and  acceleration, 
and  jerk  restricted. 

3.  The  time  of  these  cycles  should  be  minimized. 

struction  which  will  destroy  the  accuracy  of  the 
machine. 

5.  Actuator  volume  can  be  a part  of  link  structure 
volume  and  excessive  volume  reduces  dexterity  of  a 
manipualtor. 

These  requirements  reflect  the  present  state  of  the  art  in 
the  design  of  electric  motor  drives  for  robotics. 

The  study  of  actuating  systems  can  be  divided  into  two 
areas:  matching  the  load  and  actuator,  and  optimizing  the 

velocity  trajectory  using  proper  design  of  a control  system. 
The  control  theory  has  been  improved  to  upgrade  the  manipu- 
lator performance.  Luh  [131  analyzed  the  motion  control 
for  industrial  robots  to  have  minimum  positioning  error.  By 


controlling  the  robot  at  the  joint  level,  on-line  coordinate 
transformations  are  eliminated  to  reduce  the  computing  time. 
Paul  [14]  formulated  the  dynamics  equation  and  solved  the 
control  of  manipulators  so  that  the  specified  position  tra- 
jectory can  be  achieved  by  adjusting  the  feedback  gain. 

However,  the  design  of  the  electric  motor  for  the 
actuator  is  still  separated  from  the  design  of  such  system 
components  as  the  converter,  gearbox,  and  other  components 
of  the  control  system.  To  upgrade  the  whole  system,  a study 
to  match  load  and  actuator  systems  should  be  pre-executed 
e design  of  control  systems.  This  matching  load 
r system  can  maximize  the  performance  indices 
by  optimization. 


t in  Electric  M 


Although  the  design  technology  of  the  motor  has  been 

different  manufacturing  companies  because  the  design  tech- 
nology developed  by  each  company  is  kept  secret  from  the 
others.  For  example,  in  the  Gould  Electric  Company's  report 
[15],  the  computer  aided  design  procedure  of  the  induction 
motor  tells  just  the  method  of  making  the  input  list  to  the 
computer  program.  (The  major  portion  of  this  method  is 
based  on  Alger's  book  [16].)  All  the  papers  related  to  the 
design  algorithm  of  motors  published  in  the  last  20  years 
are  from  individual  scientists  rather  than  from  manufactur- 


ing companies. 


Before  explaining  each  design  method,  let  us  first 
discuss  the  general  procedure  of  optimum  design.  When  the 
motor  is  to  be  designed  or  selected,  the  type  of  motor 
application  should  be  clarified  as  show  in  Pig.  1-1.  The 
known  information  is  different  from  design  types  as  listed 
in  Table  1-1.  Normally,  the  load  duty  cycles  are  in  com- 
plicated form;  so  an  average  of  load  cycles  is  taken  to 
determine  the  rated  torque  or  power.  When  all  the  specifi- 
cations based  upon  Table  1-1  are  considered  at  the  initial 
stage  of  design,  the  dependencies  between  all  optimum  design 
methods  will  be  set  as  follows  to  organize  the  design 
structure  in  a hierarchical  manner: 

1 . The  order  of  the  design  heirarchy  is  decided  by 
the  importance  of  the  design  criteria  which  are 
stated  by  the  customer  and  by  their  involvement  in 
different  design. 

2.  The  mandatory  part  of  motor  specifications  are 
considered  first,  and  the  rest  will  be  applied  in 
the  next  stage  of  the  design. 

Based  on  the  type  of  motor  design  in  Pig.  1-1,  the 
structure  of  an  optimum  design  procedure  is  organized  as 
shown  in  Fig.  1-2.  This  figure  shows  an  example  of  TYPE  I 
motor  design.  From  the  main  part  of  input  specifications 
and  the  initial  values  of  design  variables,  the  basic  motor 
structure  is  determined.  The  design  computations  will  find 
parameters  which  are  necessary  to  compute  the  major  criteria. 
The  geometrical  constraints,  which  are  to  be  kept  strictly, 
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and  the  magnetic  constraints  will  be  checked.  If  these  con- 
straints are  satisfied,  the  design  is  compared  with  others 
by  the  major  criterion  which  is  the  cost  function  of  an 
optimum  design.  If  it  is  not  the  optimum  design,  another 
design  will  be  tried  by  adjusting  the  design  parameters  until 
the  optimum  is  found.  This  step  of  finding  the  basic  struc- 
ture is  called  the  first  step  optimization;  the  second  step 
optimization  will  find  the  detailed  design  parameters  for  the 
specific  design  and  the  performance  check. 

The  stages  development  of  design  technology  can  be 
arranged  as  follows: 

1 . Conventional  design 

2.  Initial  stage  of  computer  use 

3.  Early  stage  of  optimum  design 

4 . Optimum  design 


1.2.1  Conventional  Design 

At  the  beginning  of  this  century,  a motor  was  designed 
by  the  trial-and-error  method  [6,  16,  17];  when  the  proper 
values  of  motor  dimensions,  number  of  windings  and  the 
voltage  level  were  given,  the  performance  of  the  assembled 
motor  was  checked  and  modifications  tried  by  some  variations 
of  the  design  parameters  from  experience. 

Figure  1-3  illustrates  the  conventional  design  of  the 
DC  motor.  From  the  input  parameters — power,  voltage,  and 
speed — the  main  dimensions  of  rotor,  diameter  and  length  of 
rotor  are  determined  by  the  information  of  similar  motors. 


Then  the  magnetic  flux  density  and  magnetomotive  force  (mmf) 
are  determined  in  the  same  way.  In  armature  design,  the 
number  of  armature  windings  and  current  are  computed  from 
the  operating  voltage,  speed,  and  magnetic  flux.  After  the 
number  of  slots  per  pole  is  given,  the  armature  yoke  dimen- 
sion is  computed.  To  generate  the  required  magnetic  flux, 
the  mmf  of  field  winding  is  computed  using  the  magnetic 
circuit  theory.  Then  the  diameter  of  field  coil  and  field 
pole  dimensions  are  computed  successively.  After  all  these 
detailed  dimensions  of  motor  are  determined,  the  losses 
and  weight  of  each  part  are  obtained. 

In  this  design,  the  experience  translated  into  a number 
of  heuristic  rules  is  applied  to  make  minor  changes  which 
may  improve  the  design  or  correct  it  to  specific  require- 
ments. Thus,  a good  mathematical  model  of  the  motor  design 
is  not  required. 

However,  the  output  of  a motor  depends  on  the  dimen- 
sions of  rotor,  mmf  of  winding,  and  the  speed  of  the  motor. 
The  values  of  these  parameters  must,  therefore,  be  chosen 
to  have  optimum  cost  function,  such  as  weight,  manufacturing 
cost,  or  losses  under  a set  of  operating  constraints.  Since 
the  mmf  of  winding  can  be  expressed  as  a function  of  rotor 
dimensions  and  current  density,  these  rotor  dimensions  be- 
come the  most  important  factors  in  deciding  the  cost  func- 
tion or  performance  of  a motor. 

In  conventional  design,  these  parameters  are  determined 
by  choosing  suitable  values  obtained  from  similar  motors  or 


these  parameters 


from  experience , but 
the  optimization  which  may  b 
pose  or  cost  function.  When 
puted  from  the  information  o 
ience,  the  final  result  will 


can  be  determined  from 
different  for  a different  pur- 
all  of  these  variables  are  corn- 
similar  motors  or  from  exper- 
not  be  the  optimum  design. 


1.2,2  Initial  Stage  of  Computer  Dse 

When  the  computer  was  first  introduced,  the  majority  of 
the  earliest  problems  to  be  attempted  in  the  field  of  motor 
design  by  computer  were  of  the  formula-translation  types. 

Many  of  these  problems  and  their  solutions  have  been  described 
in  relation  to  electrical  motors.  In  these  cases,  the  com- 
puter served  as  an  aid  to  the  designer  by  allowing  a large 
number  of  calculations  to  be  made  and  more  accurate  methods 
to  be  used  and  by  disposing  of  the  routine  and  tedious 
calculations . 

Several  studies  [18-26J  are  categorized  in  this  type  of 
design.  At  this  stage,  they  translated  the  formulas  of 
motor  design  parameters  into  computer  languages.  Based  upon 
the  designer's  insight  and  previous  experience,  a set  of 
dimensions  and  parameters  can  be  chosen,  and  if  the  con- 
straints are  not  satisfied,  the  designer  must  draw  upon  his 
insight  to  modify  the  parameters  to  yield  a motor  whose 
performance  will  meet  the  specifications.  These  whole 
sequences  which  are  similar  to  those  of  conventional  design 
are  repeated  under  different  values  of  design  parameters  to 


algorithmic 


compare  the  result.  Actually,  there  were  no 
improvements  in  optimization  of  motor  design. 

1.2.3  Early  Stage  of  Optimum  Design 

As  the  computer  becomes  faster,  the  algorithm  of  find- 
ing the  minimum  point  of  cost  function  using  numerical 
methods  improves.  In  1966,  Erlicki  and  Appelbaum  [27) 
found  the  minimum  overall  cost  of  an  induction  motor  by  an 
algorithm  they  developed.  This  method  consisted  of  assuming 
numerical  values  for  the  (n  - 1)  variables  and  determining 
the  minimum  for  the  n**1  variable.  This  value  is  substituted 
in  the  equation,  and  the  minimum  is  found  for  another  vari- 
able at  (n  - 1)  constant  numerical  values  for  the  rest,  in- 
cluding the  minimum  point  found  earlier.  This  procedure  is 
repeated  until  all  n variables  are  found.  But  this  design 
problem  is  treated  as  one  with  unconstrained  minimization; 
hence,  the  results  do  not  correspond  to  reality. 

Anderson  [28]  used  Monte  Carlo  routines  to  find  the 
minimum  cost  of  a hydroelectric  generator.  After  computing 
the  material  cost  plus  the  capitalized  cost  of  losses  at  a 
certain  starting  point,  variables  on  all  design  parameters 
are  checked  to  find  the  cost  decreasing  directions.  If  the 
tries  of  some  variables  are  successful  in  reducing  cost,  the 
movement  continued  in  that  direction;  if  not  successful, 
opposite  directions  are  tried.  These  processes  will  be  con- 
tinued until  the  minimum  point  is  found.  Other  studies  (29, 


30]  are  similar  to  those  two  approaches.  The  summarized 
description  of  each  design  is  shown  in  Table  1-2. 

These  methods  were  the  first  tries  of  optimization 
before  the  nonlinear  programming  algorithm  was  widely 
developed.  Those  algorithms  took  more  computing  time 

verge  to  an  exact  optimum  point  because  of  the  discontinuity 
of  cost  functions  and  constraints.  When  the  systematic 
search  to  the  optimum  point  is  employed  by  finding  the 
effective  descent  direction  which  was  based  upon  the  non- 
linear programming  algorithm,  the  global  optimum  will  be 
achieved.  But  this  clearly  showed  the  direction  in  which 

1.2.4  Optimum  Design 

As  optimization  theory  and  the  search  algorithm  were 
more  widely  developed,  the  convergence  speed  was  increased 
by  searching  effectively  the  cost  function  through  the  de- 
creasing direction.  As  a result,  the  computing  time  is 
reduced,  and  a global  optimum  point  is  obtained. 

Schinzinger  [31]  formulated  the  transformer  design  as  a 
problem  in  nonlinear  programming.  The  cost  function  of  a 
transformer,  which  is  the  sum  of  material  cost  and  operating 
cost,  is  minimized  subject  to  constraints  such  as  efficiency 
and  temperature  rise.  He  used  the  univariate  search  which 
finds  the  descent  direction  by  changing  a single  variable  at 
a time  and  moves  to  the  descent  direction  to  further  reduce 


function . 


Ramarathnam  and  Desai  [32,  331  used  Sequential  Uncon- 
strained Minimization  Technique  (SUMT)  to  optimize  the  design 
of  the  induction  motor  with  respect  to  the  material  cost. 

They  used  external  penalty  function  to  change  a constrained 
problem  into  a unconstrained  one,  and  a first-order  gradient 
method  was  employed  to  move  into  descent  direction. 

Menzies  and  Neal  [341  formulated  the  optimization  prob- 
lem for  large  induction  motor  design  using  the  pattern  search 
and  penalty  function  method.  The  objective  function  is  the 
sum  of  material  cost  and  a specified  constrained  cost  which 
is  related  to  constraints. 

Fulton  et  al.  [351  found  the  optimum  design  in  a dif- 
ferent way  from  the  other  studies  listed  above.  They  checked 
the  variations  of  the  material  cost  of  an  induction  motor 
graphically  with  three  variables  and  found  better  design  with 
less  cost.  The  more  economical  combination  of  core  stack 
length,  wire  diameter,  and  number  of  turns  per  coil  could 
be  found  by  repeating  these  procedures.  This  method  is 
simple,  but  the  exact  optimum  value  is  hard  to  find  by  this 
method . 

Perez  and  Kassakian  [361  designed  a high-speed 
synchronous  motor  with  respect  to  minimum  size  using  the 
pattern  search  method  of  Hooke  and  Jeeves  and  an  objective 
function  which  is  a variant,  proposed  by  Menzies  and  Neal 
[341,  of  the  least  p optimization.  This  allows  the  search 


designed  by  Campbell 


al.  (37].  After  the  initial  analysis  of  finding  the  ratio 
of  inner  and  outer  diameter  of  the  motor,  the  computer-aided 
design  was  applied  to  find  the  final  optimum  design.  In 
this  analysis  step,  assumptions  were  made  in  order  to  make 
the  analytic  computation  possible , but  the  computer  could 
have  solved  the  problem  accurately  without  any  assumption. 
There  were  some  other  studies  1 38-40)  which  were  similar  to 
those  approaches  mentioned  above. 

As  can  be  seen  from  all  these  optimization  strategies 
shown  in  Table  1-2 , the  algorithms  of  finding  the  optimum 
were  improved  to  search  for  the  global  optimum  with  less 
computing  time.  The  targets  of  these  designs  were  the  gen- 
eral purpose  motors  that  satisfy  the  requirements  from  broad 
market  customers.  In  most  of  the  optimum  designs,  the  cost 
function  or  objective  function  is  formulated  using  many 
variables  and  optimized  with  respect  to  all  design  param- 

more  complicated. 

This  study  divides  the  whole  optimization  procedure 
into  several  design  stages,  and  each  design  stage  is  opti- 
mized and  linked  together  into  one  software.  The  number  of 
design  variables  is  reduced  by  expressing  some  parameters 
with  design  variables.  As  a result,  this  hierarchical  de- 
sign, with  fewer  variables,  can  noticeably  simplify  the 
design.  This  kind  of  hierarchical  design  has  the  advantage 
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of  flexibility;  a detailed  design  can  be  done  using  the  same 
main  software  program. 

This  dissertation  is  organised  as  follows.  In  Chapter 


of  motor  design  are  defined.  In  Chapter  II,  a database  is 
organized,  and  the  characteristic  parameters  of  motors  are 
compared  for  different  types  and  for  different  manufacturing 
companies.  In  Chapter  III,  the  generalized  electric  motor 
is  formulated,  and  transformer  and  synchronous  motor  designs 
are  optimized.  In  chapter  IV,  a DC  PM  motor  is  optimized, 
and  the  result  is  compared  for  two  different  materials  and 
for  two  different  structures.  In  Chapter  V,  the  design  of 
DC  PM  motor  for  multilink  manipulator  is  carried  out  to 
reduce  the  time  of  operation  using  the  constraint  of  maxi- 
mum bending  capacity.  In  Chapter  VI,  this  dissertation  is 
concluded  with  suggested  design  structure,  conclusions. 


suggestions 


future  research. 


CHAPTER  II 

ANALYSIS  OF  DATABASE  "MOTOR-USA” 


All  of  the  industrial  processes  require  the  electro- 
mechanical energy  conversion  from  a readily  available  form 
to  another  form  more  suitable  for  utilization.  There  is, 
thus,  a large  variety  of  modifications  of  electric  motors 


used  in  different  branches  of  industry.  The  real  expres- 
sions of  those  electromechanical  models  are  reflected  in 
sets  of  data  for  all  the  products  that  absorb  multi-dimen- 
sional purposes.  The  information  that  may  be  presented  for 
the  customers  is  concentrated  in  the  catalogues.  From  these 
catalogues,  one  can  extract  the  representative  character- 
istics of  those  motors.  Different  companies  submit  differ- 
ent representative  characteristics  in  their  catalogues. 
Therefore,  the  lacking  data  should  be  evaluated  analytically 
or  are  subject  to  a special  inquiry. 

An  attempt  is  undertaken  to  collect,  within  the  frame- 
work of  a single  database,  the  representative  characteris- 
tics for  a large  number  of  motors  (manufactured  in  USA) . 

The  database  is  described  in  this  study.  It  is  organized  on 


the  basis  of  the  collected  catalogues 


following 


purposes . 


manufacturing  industry  by  computing  the 


characteristic  parameters  designating  the 


lations  between  T/W  and  between  power , and 
efficiency  and  power. 

2.  To  provide  the  information  on  motor  selection  by 
comparing  the  characteristic  parameters  of  different 
types,  different  speeds  and  powers. 

3 . To  find  out  whether  the  optimism  geometry  of  the 
motor  is  actually  present  in  industry. 

This  database  is  an  integration  of  the  distinct  data  of 
all  types  of  electric  motors  and  is  stored  in  NERDC  disks. 

The  access  to  this  database  is  handled  by  SAS  software  which 
is  a computer  system  for  data  analysis  and  provides  the 
following  tools  needed  for  data  analysis  [411: 

1.  Information  storage  and  retrieval 

2.  Data  modification  and  programming 

3.  Statistical  analysis 

4.  File  handling 

So  this  database  has  the  following  advantages  among  the  gen- 
eral characteristics  of  database  [421  in  the  analysis  of 

1.  Redundancy  can  be  reduced  by  selecting  independent 
data  entries. 

2.  The  data  can  be  shared  among  different  users  by 
using  different  application  programs  in  the  SAS 

updated. 
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4.  The  number  of  words  in  the  list  of  characteristic 
parameters  can  be  easily  increased. 

2.1  Structure  of  Data 

During  the  preparation  of  the  database,  more  than  300 
manufacturing  companies  were  screened,  and  1815  motors  from 
34  manufacturing  companies  were  selected  based  upon  the 
following  selection  criteria. 

1.  The  general  purpose  motors,  except  gear  motor, 
torque  motor,  and  servo  motor,  are  chosen  as  the 
target  of  analysis. 

2.  The  information  for  weight,  power,  and  speed  must 
be  described  in  the  catalogue. 

The  information  in  this  database  includes  power,  voltage, 
current,  dimensions,  efficiency,  weight  in  AC  motors.  And 
DC  motors  add  inertia  of  motors  to  the  information  of  AC 

The  power  ranges  of  motor  application  in  common  use  are 
divided  into  five  subgroups,  as  shown  below,  because  some 
types  of  motors  have  concentrated  in  certain  power  ranges 
as  shown  in  Fig.  2-1. 


power  group 


The  RPM  groups  are  divided  into  three  subgroups  because 
in  AC  machines  the  speed  ranges  in  America  are  1200,  1800, 


respectively.  From  this  figure,  one  can  determine  the  power 
ranges  of  major  applications  in  each  type  of  motor  because 
the  number  of  catalogues  can  explain  the  number  of  motors 
manufactured  and  used  in  common  use.  For  example,  induction 
motors  are  used  in  all  power  ranges  while  synchronous  motor 

eral  purposes.  In  DC  motors,  series  and  shunt  motors  are 
being  replaced  by  permanent  magnet  motors,  and  the  develop- 
ment of  rare  earth  PM  material  will  accelerate  this  re- 
placement to  the  high  power  range  applications. 


2.2.1  The  Characteristic  Parameters 

The  performance  of  a motor  can  be  judged  and  compared 
with  other  motors  by  checking  the  specifications  and  com- 
parison criteria  as  was  shown  in  Table  1-1.  The  reference 
for  motor  evaluation  in  comparison  with  other  motors  can  be 
set  by  several  characteristic  parameters.  For  the  selection 
of  motors  in  actuators  for  robot,  aircraft,  and  automatic 


manufacturing  systems,  higher  power,  torque,  and  efficiency 
with  less  weight  are  the  primary  criteria. 

For  that  purpose,  power  per  weight  (P/W)  can  be  used, 
which  expresses  the  power  capacity  that  can  be  transformed 
from  electrical  energy  into  mechanical  energy.  In  the  same 
way,  torque  per  weight  (T/W)  presents  the  torque  capacity  of 
a motor  per  unit  weight.  This  T/W  ratio  is  similar  to  the 
output  coefficient  [6,  16,  17]  that  has  been  used  for  a long 
time  as  a first  step  in  designing  electric  motors.  The 
coefficient  is  defined  as  the  output  power  per  unit  speed 
and  unit  volume.  Since  the  weight  of  a motor  is  propor- 
tional to  its  volume,  the  output  coefficient  is  proportional 
to  its  T/W  ratio.  The  T/W  ratio  is  more  meaningful  than  P/W 
ratio  because  the  torque  of  a motor,  rather  than  its  power, 
decides  the  size  of  the  motor.  Efficiency  can  be  used  as  a 
characteristic  parameter  of  motor  because  it  is  related  to 
temperature  rise  and  rating  of  motors. 

2.2.2  Comparison  of  T/W  Ratio 

From  the  data,  the  T/W  ratio  can  be  compared  in  each 
subgroup.  In  a very  small  power  group  Cl  (0.00001-0.1  HP), 
the  bar  chart  in  Fig.  2-2a  shows  the  average  values  of  T/W 
in  each  type  of  motor  by  different  RPM  groups.  DC  PM  motor 
shows  the  highest  ratio  in  RPM  group  A because  PM  motors  can 
be  applied  at  very  low  speed.  The  induction  motor  also  shows 
higher  ratio  because  of  its  simple  structure.  When  the  motor 
is  sorted  by  torque  ranges,  the  PM  motor  shows  the  highest 


Comparison 
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other  types 


This  fact  shows  that  the  PM  motor  has  better  characteristics 
in  small  power  range  with  ferrite  and  alnico  materials , but 
when  rare  earth  PM  material  is  used,  the  difference  between 
PM  motors  and  other  types  of  motors  will  become  larger. 

Figure  2-3  illustrates  the  dependencies  between  the 
average  value  of  T/W  ratio  and  the  rated  power  for  induction 
motors  and  DC  motors.  The  plain  line  shows  the  best  values, 
and  the  dotted  line  shows  the  average  values  of  the  ratio. 

From  these  curves,  the  approximated  formulas  can  be  deter- 
mined as  shown  in  Table  2-1.  We  can  use  these  formulas  as 
references  at  the  initial  stage  of  motor  selection.  For 
example,  if  the  power  range  and  speed  are  known,  then  the 
best  and  mean  values  of  T/W  ratios  and  weights  of  motors 
manufactured  in  industry  can  be  computed  successively 
(similar  attempts  are  known  from  [7]).  The  detailed  values 
of  T/W  ratios  in  each  type,  speed,  range,  and  power  are 
given  in  Table  2.2  From  these  data,  one  can  get  more  detailed 
information  on  the  weight  and  torque  of  motor  that  will  be 


2.2,3  Comparison  by  Manufacturing  Company 

According  to  the  collected  data,  the  specifications  of 
motors  are  substantially  dispersed,  which  reflects  not  only 
the  variety  of  applicational  requirements  but  also  some 
differences  in  accepted  approaches  to  a design  procedure. 
The  trade-off  between  determination  to  provide  a number  of 


Pig.  2-3 


desirable  features  such  as  lighter  motors,  higher  efficiency, 
acceptable  manufacturing  costs,  and  customers'  technological 

To  compare  the  characteristic  parameters  of  motors  from 
different  designs,  five  manufacturers  of  induction  motors  are 
selected  in  this  analysis  with  the  following  abbreviations 
as  shown  in  the  next  figures. 

A:  Westinghouse 

C : Electrim 

D : Marathon 

The  mean  values  of  T/W  (or  P/w)  in  the  bar  chart  (Fig. 
2-4)  shows  that  motors  from  company  "B"  have  a 35%  higher 

n each  company  are  illustrated  in  Figs. 

ios  are  increased  as  the  power  becomes 
are  a couple  of  decreasing  points.  This 
does  not  mean  that  the  technology  of  motor  manufacturing 
becomes  worse  in  that  power  range.  It  does  mean  that  motor 
is  overdesigned,  e.g.,  equipped  with  a bigger  frame  and  case 
than  necessary.  This  difference  becomes  larger  in  RPM  group 
C,  which  means  that  in  high  speed,  more  different  technol- 
ogies are  applied. 


from  company  "C.”  T 
different  powers  i 
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Considering  these  factors,  the  difference  is  still 
large,  so  there  are  still  areas  to  be  improved  in  this  T/W 
ratio  which  means  the  weight  of  the  motor  can  be  further 
reduced. 

2.2.4  Characteristic  Parameters  for  DC  PM  Motor 

The  detailed  geometrical  analysis  was  applied  to  the  DC 
PM  motor  which  has  the  best  T/W  characteristics.  In  Fig. 

2-6,  the  relations  among  the  mean  values  of  T/W,  power,  and 
L/D  ratio  are  plotted.  In  each  group,  the  relation  of  T/W 

these  equations  are  drawn  in  this  figure.  As  can  be  seen 
from  this  figure,  the  ratio  group  1.3  (L/D  = 1.2-1. 4)  has 
obviously  maximum  values  of  ratio.  As  the  ratios  are 
farther  from  1.3,  the  values  of  T/W  ratios  decrease  more 
although  some  deviations  are  observed. 

The  relation  between  efficiency,  power,  and  L/D  has  a 
different  maximum  ratio.  The  efficiency  has  maximum  relation 
when  the  ratio  group  is  2.0  (L/D  = 1.8-2. 2)  as  shown  in 
Fig.  2-7  and  Table  2-4. 


2.3  Conclusions 

From  the  analyses  of  data,  the  following  conclusions 
can  be  obtained: 


applications. 


suitable 


The  best  configuration  (L/D)  that  gi\ 
performance  is  different  from  motors 
application . 


the  robotic 


The  characteristic  parameter  in  all  types  of  motors 
can  still  be  improved  by  adjusting  the  geometry  or 
by  changing  the  frames  without  introducing  special 
techniques  of  manufacturing. 


CHAPTER  III 

GENERAL  APPROACH  TO  THE  PROBLEM  OP  OPTIMUM  DESIGN 
OF  ELECTRIC  MOTOR 

3.1  Types  of  Optimum  Design 

The  optimum  design  of  an  electric  motor  will  have 
different  results  if  the  purpose  of  the  motor  is  changed. 
Table  3-1  shows  different  types  of  motor  applications. 

The  general  purpose  motor  has  continuous  ratings;  the 
output  can  be  carried  indefinitely  without  exceeding  estab- 
lished limitations.  The  optimisation  of  this  motor  is  to 
find  a motor  design  that  has  minimum  initial  cost  or  the  sum 
initial  cost  and  operating  cost.  The  former  is  dependent 
on  the  weight  of  the  motor,  and  the  latter  is  dependent  on 
the  losses  of  the  motor.  This  design  is  constrained  by  the 
performance  requirements  specified  by  National  Electrical 
Manufacturers  Association  (NEMA) . 

In  special  purpose  motors,  the  motors  for  actuators  in  a 
robotic  manipulator  [43,  44]  will  be  designed  to  have  minimum 
weight  because  a motor  with  less  weight  will  have  less  load 
inertia  reflected  to  the  next  joint  motors,  and  the  motor  will 
reduce  energy  consumption.  The  optimum  design  is  constrained 
by  higher  peak  torque  in  addition  to  the  standard  performance 
requirements.  An  electric  vehicle  [45]  is  similar  to  the 
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robotic  applications  except  that  the  power  source  is  a 
battery  and  the  operating  voltage  is  low,  which  will  make 
the  objective  function  the  sum  of  the  weight  and  the  losses 
of  the  motor.  The  crane  motor  [46)  is  another  type  of 
special  purpose  motor.  This  motor  can  usually  be  applied 
to  the  motors  which  carry  out  the  task  illustrated  in  Table 
3-1 . Crane  motor  requirements  are  often  very  severe  and 
require  a high  standard  of  reliability  to  assure  uninterrupted 
service  as  well  as  safety  to  personnel  and  equipment.  Motor 
design  must,  therefore,  have  maximum  reliability  with  in- 
creased cost  and  weight.  This  motor  will  have  an  ample  metal 
section  with  adequate  insulation  and  bracing  of  the  coil 
windings  and  connections.  The  problem  of  protecting  the  in- 
sulator can  be  achieved  by  reducing  the  losses  of  motor.  One 
can  see  how  the  difference  in  the  application  affects  the 
accepted  values  of  T/w  or  P/w. 

Ii2 General  Problem  of  Electrical  Machine  Optimization 

3.2.1  Modelling  and  Formulation  of  Magnetic  System 
The  doubly  excited  magnetic  field  system  with  two 
electrical  terminal  pairs  and  rotational  mechanical  displace- 
ments will  be  modelled  and  optimized  as  a general  form  of 
electrical  machine  as  shown  in  Fig.  3-1.  When  there  is  no 
air-gap,  this  system  will  be  changed  into  a pot  core  trans- 
former. This  system  will  be  simplified  to  have  the  following 
assumptions  without  loss  of  generality. 


*1  = f"niBrLD/2  ** 
m L1I11  + l12i2  Cos  ® 


linkage  ii 


secondary  winding 


(3.4) 


When  the  secondary  winding  is  rotated,  the  mutual  inductance 
term  is  changed  at  each  angular  position  which  generates  the 
torque  on  the  winding.  The  magnetic  flux  in  the  secondary 
winding  is  obtained  from  the  integration 

«2  =/  BrLD/2  d« 

= _fg_(Nlil  Cos  6 + "aV  <3-5> 

The  conservation  of  energy  in  the  magnetic  circuit  can 
have  the  following  relation. 


^Tldej 


e coenergy  stored  i 


The  torque  which  is  the  partial  derivative  of  coenergy  w 
respect  to  the  angular  displacement  can  be  expressed  as 


b average  torque  around  the  periphery  inside  the  stator  i 


If  this  magnetic  system  can  be  rotated  with  speed  w by 
other  additional  devices  {ex.  DC:  comutator  in  the  rotors 

AC:  rotating  magnetic  field) , the  output  power  of  this 

rotating  system  will  be 


(3.10) 


Design 


3. 2. 2.1  Input  and  design  variables 

electric  motor  was  illustrated  in  Fig.  1-2.  For  this  system 
design,  the  output  power,  operating  voltage,  and  rated  speed 
will  be  chosen  as  input  parameters,  and  the  independent 
design  parameters  must  be  chosen  carefully.  Otherwise,  more 
constraints  should  be  applied  to  this  design,  which  will  make 
the  optimum  design  still  more  complicated.  Thus,  diameter 
(D)  and  length  (L)  of  the  rotor,  window  areas  of  primary  and 
secondary  windings  (Aj , A,) , and  current  densities  of  both 
windings  (Jj,  Jj)  can  be  chosen  as  design  variables  because 
these  parameters  will  govern  the  weight,  efficiency,  and  the 
cost  of  the  motor  noticeably.  Other  parameters,  such  as 
number  of  poles,  number  of  slots  per  pole  (Nas) , and  air-gap 
length  (g)  are  assumed  to  be  fixed  for  each  motor.  And  the 
thickness  of  the  stator  can  be  expressed  as  a function  of 
design  variables. 


3. 2. 2. 2 Objective  function 

The  objective  function  of  this  system  can  be  weight,  or 

for  which  the  machine  is  going  to  be 
this  system  will  be 


used . The  weight 


Wj  = iron  weight 
= "dj { (D+2g+2tf) 2 

- (D+2g)2)L/4-Ld1(2A1+iiD2/4-2ft2)  (3.12) 

W,  = copper  weight 

= d2Aj(2L+2D+4gl  + d2A2(2L+2D)  (3.13) 

The  losses  of  this  system  will  be 

Tloss  - L1  + L2  13-14) 


= copper  losses 
= oJjAj (2L+2D+4g)  h 


(3.15) 


Cj  = core  losses  per  unit  weight 
3 .2, 2. 3 Constraints  in  optimization 

The  optimisation  of  this  system  will  be  done  under  the 
following  constraints. 

1.  Geometrical  constraints 

Available  flux  path  in  rotor  > 0 

2.  Magnetic  constraints 

Magnetic  flux  density  < saturation  density 


3.  Thermal 


Maximum  temperature  < limit  temperature 
The  geometrical  constraints  should  be  strictly  adhered  to; 
otherwise,  the  motor  cannot  be  made. 

3.2.3  Optimization  Strategy 

After  all  these  preparations,  the  objective  function 
will  be  optimized  while  maintaining  the  constraints.  The 
procedure  of  optimization  was  illustrated  in  Pig.  1-2. 

3.2.3. 1 Penalty  function  method 

There  are  many  techniques  available  for  the  unconstrained 
optimization  problems,  but  the  special  property  of  objective 
functions  of  motors  made  it  hard  to  select  the  proper  opti- 
mization algorithm.  The  characteristics  of  the  motor  objec- 
tive function  are  as  follows. 

1 . The  number  of  turns  of  the  winding  have  to  be  an 

2.  The  variables  of  all  equations  are  linked  together, 
so  it  is  very  difficult  to  differentiate. 

3.  The  variables  are  to  be  chosen  under  the  constraints. 

4.  As  a result,  the  objective  function  is  not  oontin- 

Therefore,  the  algorithms  using  derivatives  [47-491  are  not 
suitable  for  this  purpose.  A penalty  function  method  that 
transforms  the  constrained  problem  into  a sequence  of  uncon- 

There  are  several  reasons  for  the  appeal  of  penalty  function 


Fig.  3-2,  is  that  the  s 


Pig. 


Illustration  o£  penalty  function 


Since  this  equation  does  not  allow  any  constraint  to  be  vio- 
lated, it  requires  a feasible  starting  point  for  the  search 
toward  the  optimum  point.  However,  in  many  engineering 
problems,  it  may  not  be  so  difficult  to  find  a point  satis- 
fying all  the  inequality  constraints.  If  0(x,rk)  is  mini- 
mized for  a decreasing  sequence  of  values  r , the  uncon- 
strained minima  will  converge  to  the  solution  of  the  original 
problem  stated  in  eq.  (3.17)  as  shown  in  Fig.  3-2. 

3. 2. 3. 2 Search  algorithm 

To  achieve  a rapid  convergence  to  a global  minimum  of 
a N variable  function,  the  pattern  search  method  is  used  in 
this  study.  The  pattern  search  of  Hooke  and  Jeeves  (51)  is 
a sequential  technique,  each  step  of  which  consists  of  two 
kinds  of  moves — an  exploratory  move  and  a pattern  move.  The 
exploratory  move  starts  X°  as  the  initial  point,  and  all 
coordinate  directions  are  searched  one  at  a time  to  find  the 
descent  direction,  if  a positive  displacement  of  any  point 
does  not  yield  a reduction  of  the  objective  function,  the 
opposite  direction  is  tried.  At  the  end  of  a series  of 
exploratory  moves,  a new  point  X1  can  be  determined.  Now  a 
pattern  search  is  made  along  (X*  - X^)  direction  and  is  con- 
tinued until  a minimum  is  found  in  that  direction  as  shown 
in  Fig.  3-3.  At  the  end  of  this  pattern  search,  the  pro- 
cedure is  restarted  until  a global  minimum  is  found. 

For  the  minimization  technique  in  one  direction, 
Fibonacci  search  is  used.  Fibonacci  search  [521  reduces 


the  interval  o£  uncertainty  about  where  the  minimum  lies  by 
eliminating  the  appropriate  interval  after  comparison  of  two 
function  values  in  two  Fibonacci  points.  The  detailed  flow 
chart  of  computer  program  will  be  shown  in  Appendix  C.  The 
examples  of  optimum  design  will  be  presented  in  the  next 
sections . 


3.3  Transformer  Design 


cussed  here 


procedure  and  numerical  result  w 


id  the  detailed  derivations  will  be  in  Appendix 
shows  the  geometry  of  c-core  transformer 
to  be  optimized. 


3.3.1  Preparation  for  Optimization  [111 

The  input  parameters,  design  variables,  objective  func- 
tion, and  constraints  are  expressed  as  in  Table  3-2.  All 
these  parameters  are  expressed  as  functions  of  design 
variables. 


3.3.2  Optimization 

The  optimum  design  procedure  of  the  transformer  is 
similar  to  Fig.  1-2  except  that  the  input  and  design  param- 
eters are  different.  This  design  will  determine  the  main 
design  parameters,  such  as  the  geometry  (Lc,  Dm,  Dc) , window 
area  and  current  density.  The  detailed  design  can  find  the 
diameter  of  winding  and  number  of  turns.  For  the  optimiza- 
tion of  the  transformer,  pattern  search  was  used. 


Table  3-2  Transformer  design 


Parameters 

Descriptions 

Input  parameters 

1. Induced  output  VA  : S 
2. Operating  voltage  : V 

Design  variables 

1.  Length  of  winding  coil  : L 

2.  Width  of  square  coil  : Dm 

4.Current  density  : J 

Weight  objective 
function 

Wcp^eight  of  copper  winding 
Wir=weight  of  core 

Loss  objective 

Lcp=copper  winding  losses 

PLo(Dc+Dm)  ( D^- Dn>)  ^ ' "^KuK?2^  ’ 

= (L  D ) 0 ' 25 

-L • w . {LJ  :core  losses  per 

lr  lr  unit  weight) 

Constraints 

1. Operating  magnetic  flux  density 
2. Maximum  temperature  rise  of  winding 

The  weight  function  and  losses  function  of  the  trans- 
former are  plotted  in  three  dimensions  as  in  Figs.  3-5  and 
3-6.  These  plots  show  that  there  exists  an  obvious  global 
optimum  point  for  each  objective  function.  The  numerical 
results  in  Tables  3-3  and  3-4  show  that  the  shape  of  the 
transformer  in  optimum  design  is  not  changed  although  the 
output  VA  is  increased. 

3.4  Synchronous  Motor  with  PM  Rotor 

3.4.1  Preparation  for  Optimization  [111 

The  top  view  of  a synchronous  motor  with  a PH  rotor  is 
shown  in  Fig.  3-7.  The  input  parameters,  design  variables, 
objective  function,  and  constraints  in  the  optimum  design  of 
a synchronous  motor  with  a PM  rotor  are  shown  in  Table  3-5 
and  the  detailed  derivations  are  in  Appendix  B. 

3.4.2  Optimization 

The  optimum  design  procedure  was  shown  in  Fig.  1-2,  and 
this  design  will  determine  the  main  design  parameters  such 
as  diameter  and  length  of  rotor,  window  area  of  armature 
winding  and  current  density  of  armature  winding.  The  opti- 
mization was  done  using  pattern  search  as  in  transformer 

The  numerical  result  of  optimum  design  for  minimum 
weight  is  presented  in  Table  3-6.  The  same  result  with 


of  a transformer 


Pig. 


Weight  function 


synchronous 


Parameters 

Description 

Input  parameters 

1. Output  power  : PQU^ 

2. Operating  voltage  : V 

4. Number  of  poles  : Nq 

Design  variables 

1.  Diameter  of  rotor  : D 

2.  Length  of  rotor  : L 

3.  Armature  window  area  : A 

4.  Current  density  : J 

Weight  objective 

=DacL  [(it  (D+2g+2t)  2-it  (D+2g)  2}/4 

Waw=armature  winding  weight 
_6mDcpnasNapAaLks 

Constraints 

1.  Width  of  armature  teeth  : Wa(;>  0 

2.  Magnetic  flux  density 

3.  Window  area  of  armature  winding 

4. Full  load  temperature  rise 5 r 
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transformer 


changed;  in  other  words,  the  ratio  between  diameter  and 
length  of  rotor  is  the  same  for  a series  of  different  power 


3.5  Conclusions 

Prom  all  these  analyses  and  results  of  optimization, 
we  can  get  the  following  conclusions: 

1.  All  machines  have  their  own  optimum  configurations. 

2.  The  optimum  configuration  is  not  changed  although 
the  power  rating  is  changed. 


CHAPTER  IV 

STRUCTURE  OF  DESIGN  PROVIDING  MINIMUM 
WEIGHT/LOSSES  OF  ENERGY 

The  general  problem  of  electrical  machine  optimization 
was  shown,  and  the  simplified  models  in  transformer  and 
synchronous  motor  were  optimized  in  Chapter  III.  This  chap- 
ter will  optimize  the  DC  PM  motor  with  respect  to  weight, 
losses  of  energy  and  combination  of  both  objective  functions 
in  the  normal  structure  of  a motor.  The  motor  with  an 
inside-out"  structure  will  be  optimized  to  compare  the 
advantages  of  this  structure. 

4.1  The  Specific  Characteristics  of  DC  PM  Motors 

Conventionally  wound  field  DC  motors  have  been  the  main 
DC  energy  conversion  devices  since  the  discovery  of  electro- 
mechanical motion.  Recently,  the  permanent  magnet  (PM) 
motors  challenged  the  traditional  wound  field  DC  motors  in 
all  sizes  up  to  5 HP  (53) . For  the  last  thirty  years,  the 
use  of  PM  fields  in  small  DC  motors  has  been  increasing  to 
the  point  where  a majority  of  the  motors  are  now  changed 
into  PM  motors.  Furthermore,  the  development  of  the  new  PM 
materials  also  evolved  and  broadened  their  applications. 


Prior  to  the  PM  motors,  the  three  popular  types  of  DC 
motors  for  industrial  applications  were  shunt,  series,  and 
compound  motors.  Designers  often  made  compromises  in  their 
equipment  in  order  to  use  one  of  these  three  motors.  For 
example,  gears  and  pulleys  were  often  added  to  reduce  the 
speed  of  a series  motor  enough  to  make  the  portable  equip- 
ment practical. 

For  the  proper  selection  or  design,  it  is  necessary  to 
understand  the  characteristics,  advantages,  and  shortcomings 
of  each  type  of  motor.  Table  4-1  shows  the  collected  infor- 
mation of  the  characteristics  of  wound  field  motors.  As 

a number  of  shortcomings. 

To  see  the  superior  properties  of  PM  materials,  let  us 
make  a comparison  between  the  PM  and  wound  field  methods  as 
illustrated  in  Fig.  4-1  and  Table  4-2. 

The  volume  of  field  part  is  proportional  to  the  dimen- 
sion, D,  in  wound  field  method,  and  is  proportional  to  the 
square  of  core  dimension,  D2,  in  PM  field  method.  The 
coefficient  of  volume  equation  in  the  PM  field  is  much 
smaller  than  that  in  the  wound  field  method.  Thus,  the 
volume  of  the  PM  motor  is  smaller  than  that  of  the  wound 
field  motor  in  small  power  applications.  If  the  power  is 
very  low,  the  volume  of  the  PM  motor  is  much  smaller  than 
that  of  wound  field  motor.  This  coincides  with  the  fact 


applications 


confined 


Comparison  of  wound 


Kinds 

Advantages 

Disadvantages 

“ 

1. Relatively  good  speed 
regulation  up  to  250  % 
•Variety  of  control 
•Capability  of  high 

1. Small  starting  torque 
(armature  reaction) 

motor3 

L. Highest  starting  torqu 
!. Higher  power  per  weigh 

3. Wide  speed  capability 

2. Over-speed  at  light 

Compound 

•Better  than  pure  shunt 

1 -Complicated  control 

Fig. 


field  methods 


Table  4-2  Comparison  of  PM  and  wound  field  methods 


to  smaller  power  ranges:  the  motors  £or  robotic  applications 

are  within  these  power  ranges. 

We  can  also  see  from  this  table  that  as  the  magnetic 
flux  density  increases,  the  volume  of  the  wound  field  motor 
is  increased  quite  linearly,  but  in  PM  field  motors,  the 
volume  is  increased  much  faster.  So  the  operating  point  of 

product  point. 

The  advantages  of  the  DC  PM  motor  over  the  wound  field 
motor  are  listed  below  153-55]  with  the  result  of  comparison 
from  the  database  as  shown  in  Table  4-3: 


for  the  smaller  outside  diameter  are  that  there  is 
no  bulky  field,  and  the  magnets  can  be  assembled 
by  cementing  the  segments  to  the  shell. 

2.  Weight:  A weight  reduction  up  to  30%  can  be  ex- 
pected in  the  motor  due  to  the  smaller  diameter 

base,  the  weight  of  the  PM  motor  is  46%  less  than 

3.  Efficiency:  An  increase  in  the  efficiency  by  up 

to  15%  occurs  because  the  field  is  generated  by  a 
permanent  magnet  and  does  not  have  to  be  energized. 


Table 


Ef  f iciency (% ) 
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4.  Temperature:  Higher  efficiency  means  cooler 

operation  for  the  same  power,  so  the  PM  motor  con- 
cept is  suitable  for  the  equipment  where  heat 
losses  are  problems  in  operation. 

5.  Armature  reaction:  The  PM  pole  piece  has  a per- 

meability very  close  to  that  of  air.  And  when  the 
armature  reaction  flux  leaves  the  armature  in  a PM 
motor,  the  flux  lines  are  resisted  by  a low  per- 
meability material.  In  other  words,  the  PM  motor 
has  a very  large  effective  air-gap.  So  the  armature 
reaction  does  not  distort  the  main  field  at  the 
overload  condition  in  a PM  motor  as  it  does  in  a 


wound  field  motor.  As  a result,  i 
tively  high  starting  torque  (650% 
and  the  performance  of  a PM  motor 
over  the  operating  range. 


::  rated  torque) , 

^ a straight  line 


Considering  all  these 


conversion  t 


4.2  Comparison  of  Permanent  Magnet  Materials  (56-591 

The  permanent  magnet  materials  have  been  developed  to 
have  higher  remanent  flux  density,  higher  coercive  force,  and 
higher  energy  product.  There  are  three  major  PM  families  in 
use— alnico,  ferrites,  and  rare  earth  magnet.  We  will  con- 
sider the  characteristics  of  each  PM  material. 


Demagnetization  characteristics  o£  those  PM  materials 
most  commonly  used  in  rotating  electrical  machines  are  shown 
in  Pig.  4-2  together  with  some  new  materials  of  interest. 
Table  4-4  illustrates  the  comparison  of  PM  materials. 

Alnico  magnets  have  been  traditionally  employed  in 
electrical  machines.  The  high  remanence  type,  Alnico  5-7, 
offers  potentially  higher  flux  density  with  good  tempera- 
ture characteristics.  However,  the  very  low  coercive  force 
requires  a very  long  magnet.  Furthermore,  the  knee  in  the 
curve  occurs  only  slightly  below  the  remanent  flux  density. 
This  will  result  in  the  operating  point's  being  forced  down 
the  steep  part  of  the  demagnetization  curve  with  only  a small 
additional  demagnetizing  field  due  to  an  armature  reaction. 
Another  problem  concerning  this  low  coercive  force  is  that 
this  magnet  must  be  remagnetized  after  assembling  it  to  the 


Ceramic  ferrites  are  most 

and  ease  of  fabrication  into  tt 

the  very  low  permeability  ensui 

as  with  Alnico' s.  Ferrites  mas 
bling  because  of  high  coercive 
density  is  very  low,  and  they  , 


commonly  used  in  the  produc- 
cost  advantage,  availability 
le  stator.  The  demagnetiza- 
: most  of  the  quadrant,  and 
:es  that  recoil  property  will 

r be  magnetized  before  assem- 


i frequently  designed  t 


operate  w 


ve  the  maximum  energy  point, 


Fig. 


Demagnetization  characteristics  of  PM 


developed  about 


Rate  earth  PH  materials  [60-61],  were 
20  years  ago,  and  their  properties,  which  are  almost  the 
best  in  all  areas,  are  still  improving.  The  properties  of 
these  materials  which  intrigue  the  motor  designer  most  are 
the  following. 

1.  Highest  energy  product:  144  kj/m 

2.  Highest  coercive  force:  388  kA/m 

3.  Highest  remanent  flux  density:  0.85  T 

4.  High  curie  temperature:  775“C 

The  only  disadvantage  is  the  price  of  materials,  which  limits 
the  design  criterion  of  this  type  to  operation  at  the  maxi- 
mum energy  point  to  minimize  the  magnet  volume.  With  this 
rare  earth  material,  it  is  now  possible  to  place  the  PM  field 


Figure  4-3  shows  the  normal  structure  and  "inside-out" 
structure  of  the  DC  PH  motor.  The  “inside-out"  structure  is 
the  inverted  structure  where  the  armature  part  is  in  the 
stationary  part  and  PH  is  in  the  rotor  as  shown  in  the  fig- 
ure. The  difference  in  detailed  design  will  be  discussed 

Based  upon  the  general  design  structure  in  Fig.  1-2, 
the  optimum  design  procedure  is  organized  as  shown  in  Fig. 
4-4  which  illustrates  the  optimization  algorithm,  SUMT  and 
pattern  search,  discussed  in  Chapter  III.  Most  of  the 


Pig.  4-3  Normal  structure  and  "inside-out" 
structure  of  DC  PM  motor 


Pig.  4-4  Basic  structure  of  optimum  motor  design 


procedures  are  the  same  as  the  general  structure  shown  in 
Fig.  1-2.  This  figure  explains  the  first  step  design  which 
decides  the  optimum  motor  configuration. 

After  the  input  parameters  of  DC  PM  motor  design  are 
given,  which  are  the  specifications  of  motor  performance, 
the  magnetization  and  core  loss  curves  for  selected  core  are 
determined  using  the  Chebyshev  polynomial  approximation  [62] . 
When  the  SOMT  is  started  under  given  penalty  parameter  r , 
the  pattern  search  which  consists  of  an  exploratory  move  and 
a Fibonacci  move  was  used  to  find  the  minimum  weight  or  loss 
objective  function.  After  the  series  of  minimization  pro- 
cedures, the  convergence  criterion  is  checked.  If  the 
criterion  is  satisfied,  the  optimum  design  parameters  can  be 
obtained;  otherwise,  this  optimization  procedure  (#1  in 
Fig.  4-4)  will  be  reiterated  until  the  criterion  is  satisfied. 

With  these  obtained  design  parameters,  the  temperature 
rise  of  armature  winding  is  to  be  checked  because  the 


armature  winding  is  the  source  of  heat  and  has  the  highest 
temperature  in  the  motor.  If  the  temperature  rise  is  not 
the  same  as  the  maximum  limit,  the  current  density  is  ad- 


justed and 


the  optimization  procedure  (42  in  Fig.  4-4)  will 
again  until  it  reaches  the  temperature  limit, 
will  consider  the  actual  design  of  the  DC  PM 
normal  and  "inside-out”  structure.  Although  both 
are  different  in  configurations,  their  basic 


magnetic  circuits 


Pig. 


Normal  Structure 


4 -3. 1.1  Selection  of  design  parameters 

The  flow  chart  of  computing  the  objective  function  which 
illustrates  the  design  procedure  is  shown  in  Fig.  4-6.  The 
input  parameters  are  power  (PQut) , voltage  (V) , rated  speed 
(wj , and  number  of  Poles  (N  ) for  general  purpose  motors. 

The  motor  for  the  robotic  actuator  will  have  rated  torque 
(Tr)  and  number  of  poles  (N  ) as  input  parameters. 

It  is  very  important  to  select  the  design  variables 
properly,-  otherwise,  the  optimization  problem  becomes  compli- 

design  parameters  finally  selected  in  DC  PM  motors  are 
shown  below: 

1.  Rotor  diameter:  D 

2.  Rotor  length:  L 

4.  Width  of  magnet: 

5.  Thickness  of  field  yoke:  t^ 

6 . Current  density:  J 


4. 3. 1.2  Induced  voltage  computation 

From  the  given  voltage  and  power,  the  diameter  of 
winding  can  be  determined  using  the  following  basic  equa- 


(4.2) 


winding 


Ja  = current  density  in  armature  winding 
a = number  of  parallel  paths 


0 = resistivity  of  winding  material 
Mjt  = mean  turn  length  of  armature  winding 
e window  area 
of  armature  slots  per  pole 
f pole  pairs 


Substituting  eg.  (4.4)  into  eg  (4.2)  yields 


copper  losses 


From  eg.  (4.5),  the  diameter  of  winding  can  be  computed  as 


After  the  diameter  of  armature  winding  is  obtained,  the 
induced  voltage  and  magnetic  flux  can  be  computed.,  respec- 
tively. 

Eg  = V - IaRa  (4.7) 

* = 55~w  (4'81 


= mhas  Integer  (Aa/ (nD^/4) ) (4.9) 

- angular  speed 


e yoke,  and  field  yoke  can  be  expressed  as 
f geometrical  parameters  respectively. 


Rm  ’ uruOAm 


1.11) 


(4.13) 


U3u0”fy 
= effective  a 


= 9Da3  + w„/(l  + 5g) 

' V Aay’  At-  Afy  " cross  sectional  area  of  air-g 
magnet,  armature  yoke,  teeth, 

a thickness  of  magnet 

U2U3  = permeability  of  teeth,  armature  yoke,  and 


e magnetic  circuit  equation  c 
= *(2R  + 2R  + 2Rf)  + <R 


a expressed  a 


Prom  eq.  (4.10)  - eq.  (4.16),  the  thickness  of  permanent 
magnet  that  can  generate  the  magnetic  flux  will  be  deter- 


(4.17) 


* (2R 


g + 2Rt  = Ray/2) 

■ ♦(ir/(2u,p0  A.  ) + 2/(ii  ii 


4 . 3 . 1 . 4 Approximations  o 

To  determine  exactly  the  core  losses  an 
ration  curve,  the  Chebyshev  polynomial  and  c 
are  used.  The  results  of  approximations  are 
4-7  and  4-8.  The  actual  data  represented  b 
close  to  the  approximated  data. 


L.S  Objective  function 

The  objective  function  o 
use  the  weight,  losses. 

Armature  part 


' combination  a 


le  raagneti- 
3wn  in  Figs. 


Wflw  = armature  winding  weight  = DCpMjtNAa 
Wai  - insulator  weight 


- yoke  weight  = 


W£p  = PM  weight  = NpDpmvpm 
W£y  = yoke  weight 

■ -i-j — ( (D+  2g+2Tm  + 2t£)2  - (D+2g+2T  )2) 

Wt  - Wa  + W£  (4.20) 

2 . Losses 


(4.22) 


Lafc  = teeth  core  losses  - Lg£Wa£ 
Lay  = yoke  core  losses  - Ls2way 


Total  losses 


Tloss  = Lcp  + Lcr  «'23> 

can  be  expressed  using  weighing  factor  A. 


4. 3. 1.6  Constraints 

The  constraints  for  the  optimum  design  of  a DC  PM 
motor  which  are  similar  to  those  in  Chapter  III,  are  listed 

1.  The  width  of  armature  teeth  and  armature  window 
area  are  positive. 


e magnetic  flux  density  i 
d field  yoke  is  less 


Among  those  constraints,  the  temperature  rise  is 
easily  because  the  heat  dissipation  is  related  t 
structure,  environmental  conditions,  and  time, 
be  discussed  in  detail  in  the  next  section. 


4 . 3 . 1 ■ 7 Thermal  computation 

1.  Heat  transfer  between  rotor  and  stator.  The  heat 
transfer  characteristics  of  air-gap  in  electrical  motors  ar 


the  least  studied  of  the  many  heat  flow  paths  in  an  elec- 
trical rotating  machine.  Martinelli  (64)  has  derived  an 
expression  for  the  heat  transfer  to  a fluid  flowing  between 
two  parallel  plates  with  a constant  heat  rate  at  each  plate. 
But  that  expression  is  valid  only  when  the  fluid  is  being 
symmetrically  heated  from  each  plate.  In  our  case,  the 
fluid  is  heated  from  one  side  and  cooled  from  the  other. 

of  the  Nusselt  number  for  heat  transfer  across  air-gap  and 
compared  this  with  other  experimental  data.  His  modified 
formula  will  be  used  in  this  study. 

Now  consider  the  concentric  cylinder  flow  in  which  no 
axial  flow  occurs  as  shown  in  Fig.  4-9,  which  illustrates 
the  basic  structure  of  the  motor  for  this  thermal  analysis. 
In  this  case,  the  flow  is  actuated  only  by  the  rotation  of 
one  of  the  cylinders  and  can  be  characterized  by  a Reynolds 
number  and  by  a dimensionless  curvature  factor.  The 


R _ (equivalent  diameter) (mean  velocity) 


(4.25) 


u - peripheral  velocity  and  rotor 
The  dimensionless  Nusselt  number  for  heat  transfer  across 


(65) 


C£  : Field  part  heat  capacity 
Ca  : Armature  heat  capacity 


(4.26) 


10[Pr+ln(l+5Pr)  + 0.51nHeg(cf/2)°-5/120] 


f = friction  factor 


Reg  - Reynolds  number 

cf  = skin  friction  coefficient 

The  relation  between  skin  friction  coefficient  and  the 
Reynolds  number  approximated  from  the  measured  data  is  to 


The  heat  transfer  coefficient  by  convection  across  an 


e power  transfer  by  convection  a 


transfer 


by  radiation 


transfer 


2.  Heat  transfer  from  motor  to  ambient  air  [67] . For 
convection  in  still  air,  the  heat  transfer  coefficient  can 
be  used  from  [681 , converted  to  the  units  used  in  this 


= equivalent  diameter  of  : 
- air  film  temperature 


t dissipation  from  t 


And  the  heat  dissipation  by  radiation  is 

The  total  heat  transfer  from  the  surface  of  the  motor  to 
the  environmental  air  is  then 


(4.34) 


3.  Temperature  rise  computation.  When  all  these  heat 
dissipations  are  computed,  the  temperature  rise  of  armature 
winding,  Tfl,  and  field  part,  T-,  can  be  determined  by 
solving  the  thermal  differential  equations. 
dT 

“dt  = ,P1  ‘ Pa,/Ca  (4.3S) 

dTf 

-g§  - CPa  - Pc)/Cf  (4.36) 

These  equations  are  the  nonlinear  differential  equations; 
so  the  numerical  method,  Runge-Kutta  predictor  corrector 
method  [691  can  solve  these  equations.  The  motors  for 
robotic  applications  have  a short  intermittent  duty  cycle 
(Chapter  1X1);  so  the  temperature  rise  is  computed  on  a 
one-hour  time  period,  which  is  different  from  the  continuous 
rating  in  general  purpose  motors  shown  in  Table  3-1. 

4.3.2  Inside-out  Structure 

Since  the  basic  structures  of  normal  and  "inside-out" 
are  quite  similar,  this  section  will  discuss  just  the 
differences  between  these  two  structures. 

The  thickness  of  the  PM  magnet  in  this  structure  of 

tion  except  for  the  difference  of  reluctances  by  the  change 
of  structure.  The  changed  reluctances  are  armature  and 


field  yoke  reluctances: 


(4.37) 


= n{D  + 2(g  + D ) + T ) 


= thickness  o 


Using  eq.  (4.37)  and  eq.  (4.38),  the  thickness  o 
can  be  expressed  by  the  following  equation: 


e weight  objective  function  o 


way  * — ( (Dt2g+20as+Ta)  2 - (D+2g+2Das) 2) 

Wat  = — ( (0+2g+2Das)  2 - (D+2g)2)-N  NasDjir(5asMasI" 


The  losses  of  this  motor  are 


s different 
eat  generation 


The  heat  dissipation  of  this 
from  that  of  a normal  structure  m 
in  the  armature  is  dissipated  to  the  field  part  and  air  with 
a larger  surface  area  so  the  effective  heat  dissipation  area 
is  greatly  increased  (Fig.  4-10) . As  a result,  the  thermal 


4.4  Optimization 


the  importance  of  weight  in  manipulator 
in  Chapter  I . Based  on  the  design  pro- 
in Fig.  4-4,  the  computer  program  is 


Fig. 


organized  and  the  weight  objective  function 
using  SUMT  and  pattern  search  method.  This 
will  be  explained  using  the  structure  of  a 


is  optimized 
design  procedure 
computer  program. 


4.4.1  Computer  Program 


Main  program 
main  computer  program  v 
Determine  the  purpose  c 
is  for  general  purpose 


the  following  tasks: 
multilink  manipu- 


Read  in  the  input  data  and  initial  feasible  values 
of  design  variables. 

Find  the  magnetization  and  core  loss  curves  of 
selected  core  material. 

Control  the  SUMT  by  the  penalty  factor. 

Check  the  temperature  rise,  adjust  the  current 
density,  and  repeat  above  procedures  3 and  4 until 
the  temperature  rise  is  the  same  as  the  limit. 

Print  the  specifications  c 

’ Subprogram  "PATTRN" 


e designed  n 


The  subprogram  "PATTRN”  will  find  the  minimum  point  by 
the  following  procedures: 

1.  The  subprogram  "explor"  will  find  the  descent 

direction  by  trying  all  the  directions  of  design 
variables  that  can  reduce  the  objective  function. 
If  the  search  could  not  reduce  the  objective 


function,  the  new  values  of  design  are  replaced 
by  the  former  values  before  the  search. 

2.  If  the  exploratory  move  fails,  the  search  step 
distances  will  be  reduced  and  the  exploratory  move 
tried  until  the  descent  direction  is  found. 

3.  The  pattern  move  will  be  done  by  the  subprogram 
"FIBONA."  This  subroutine  will  find  the  minimum 
point  using  a Fibonacci  search  which  is  an  effec- 
tive line  search  method. 

4 . The  convergence  criterion  can  be  satisfied  if  the 
three  successive  results  of  pattern  search  are 
within  a satisfactory  limit,  otherwise  the  above 
procedure  will  be  reiterated  until  the  criterion 
is  satisfied. 


Subprogram  " 


The  objective  function  can  be  determined  from  this 
subprogram  "DESIGN”  which  will  find  the  performance  param- 
eters and  check  the  constraints.  The  computing  procedures 
are  shown  in  Fig.  4-6.  The  detailed  flow  charts  of  optimiza- 
tion and  computer  program  are  in  Appendices  C and  D. 


4.5  Result  of  optimization 
4.5.1  Convergence  Check 

The  penalty  function  method  and  pattern  search  method 
used  in  this  computer  program  of  optimum  design  show  rapid 
convergence  to  the  optimum  values.  Figure  4-11  shows  the 
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convergence  speed  of  optimum  design  of  a 5 N-M  motor . After 
six  iterations,  all  of  the  design  parameters  are  very  close 

tinued  for  up  to  15  iterations  to  increase  the  accuracy  of 
the  results.  Optimization  with  different  starting  points 
was  tried,  and  the  results  are  within  0.2%  difference. 


5 . 2 Computer  Output 

The  specifications  of  the  designed  m 
th  the  format  shown  below: 

Result  of  optimization 

Torque  (N-m)  

Number  of  poles  

Power  (HP)  

2.  Design  Parameters 

Rotor  diameter  (cm)  

Rotor  length  (cm)  

Magnet  width  ( cm)  

Stator  thickness  (cm)  

Current  density  (A/cnr)  . . 

3 . Electrical  Specifications 

Number  of  slots  per  pole  . 


Total  number 


Resistance  of  armature  coil  (H) . 
Inductance  of  armature  coil  (mH) 

Induced  voltage  (volt)  

Magnetic  Specifications 

Flux  (WB)  

Flux  density  (TESLA) 

Magnet  


Field  yoke  

Loss  (watt) 

Copper  loss  

Mechanical  loss  . . . . 

Total  loss  

Efficiency  (»)  

Armature  (ROTOR)  ... 
Winding  


Field  (STATOR) 


Magnet 
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Total  weight  

7.  Acceleration  (rad/sec') 

inertia  (kg-m2)  

8.  Performance  Constant 

Voltage  constant  (Ke:V. sec/rad) 

Torque  constant  (Kf:N-m/A)  

Electrical  time  constant  (msec) 

Mechanical  time  constant  (msec) 

4.5.3  Numerical  Results  with  Different  Torqui 

Table  4-5  shows  the  numerical  results  of  optimum 
design  of  weight  minimization  (X  = 0)  in  a series  of  dif- 
ferent torques.  In  lower  torque  ranges,  the  motor  is  very 
small,  and  it  is  hard  to  satisfy  the  following  constraints 
with  the  small  rotor  diameter: 

1.  The  width  of  slot  in  the  armature  is  positive  to 
put  the  armature  winding  in  it. 

2.  The  tooth  width  in  the  armature  is  positive  and  is 
larger  than  some  value  to  make  the  magnetic  flux 
density  in  the  tooth  less  than  saturation  density. 

As  a result,  the  diameter  of  a rotor  is  to  be  increased  to 
satisfy  those  constraints.  As  the  torque  becomes  larger, 
those  constraints  can  be  satisfied  easily  with  this  diam- 
eter. When  the  torque  of  the  motor  is  larger  than  5 N-m 


4-5.4  Comparison  with  Different  Material  and  Structure 

The  DC  PM  motor  is  designed  in  different  magnet  mater- 
ials, ferrite  and  SmCo_ , and  in  different  structures,  normal 
and  “inside-out.”  The  results  of  these  designs  are  illus- 
trated in  Table  4-6.  Prom  this  table,  one  can  see  the  follow- 
ing differences: 

1.  In  normal  structure,  the  SmCo.  motor  shows  30%  less 
weight,  3%  higher  efficiency,  and  35%  higher  no 
load  acceleration  than  does  the  ferrite  motor. 


4 . 6 Conclusions 

Based  upon  the  analysis  and  results  of  DC  PM  motor 
optimum  design,  we  can  get  the  following  conclusions: 

1 . The  DC  PM  motor  has  better  characteristics  than 
the  wound  motor,  as  shown  in  Table  4-6. 

2.  The  DC  PM  motor  has  a unique  shape,  but  the 
constraints  modify  the  shape  to  make  a feasible 

3.  The  most  desirable  material  for  an  actuator  in 
robotic  application  is  a rare  earth  PM  material. 

4 . The  suggested  structure  of  a motor  for  robotic 
application  is  the  "inside-out"  structure  with  a 
rare  earth  PM  magnet. 


of  optimum  design  is  developed 


based  upon  the  requirements  of  a TYPE  I design 


CHAPTER  V 
DESIGN  OF  DC  PM 
HE  ENHANCED  TIME 


' OPERATION 


Previously,  the 
presented  applicable 
to  be  equipped  by  a n 


electric  motor  design  was 
ise  when  the  working  machine 
already  designed  or  existed. 


Now  we  will  consider  an  important  case  when  the  design  of 

Based  on  the  results  of  optimum  design  of  a DC  PM  motor 
in  Chapter  IV,  the  structure  of  design  will  be  modified  to 
satisfy  the  requirements  of  enhanced  time  of  operation.  The 
time  of  operation  can  be  reduced  by  two  different  ways: 
matching  the  load  and  actuator,  and  optimizing  the  velocity 
trajectory  via  proper  design  of  a control  system.  This  study 
will  be  focused  on  the  time  reduction  by  matching  the  load 
and  actuator  conditions  presuming  a possibility  to  provide 
a simultaneous  design  of  mechanical  and  electrical  parts. 

We  assume  that  the  minimum  time  velocity  trajectory  is 
obtained  (which  is  known  to  be  a triangle  (701),  and  the 
load  parameters  such  as  load  inertia,  Jj,  angular  displace- 
ment, 0, , are  given  as  input  information.  The  design  vari- 
ables of  this  system  will  be  motor  torque,  T , and  gear 
ratio.  Under  given  torque,  the  inertia  is  determined  from 
the  optimum  design  of  motor  for  minimum  weight.  The 


optimization  of  finding  motor  torque  and  gear  ratio  will 
reduce  the  operation  time.  If  the  weight  of  motor  is  to  be 
minimized  under  each  different  value  of  motor  torque,  then 
this  time  reduction  will  require  a big  computer  program  and 
will  take  much  computing  time.  However,  if  the  motor  inertia 
can  be  expressed  as  a function  of  motor  torque  from  the 
result  of  weight  minimization,  then  this  problem  will  become 
much  simpler.  Thus,  the  hierarchical  design  structure  of 
time  minimization  in  actuator  system  is  organized  as  shown 
in  Fig.  5-1. 


en  the  velocity  trajectory  and  input  parameters  a 
the  time  of  operation  can  be  reduced  as  follows: 
Find  the  functional  relationship  between  the  mo 
inertia  and  torque. 


reduce  the  operation  time  ar 
maximum  bending  capacity  of 


3 . Compute  the  maximum  speed  and  maximum  power  in 
that  given  velocity  trajectory. 

4 . Since  the  operating  voltage  is  a function  of  mag- 
netic flux  that  is  determined  in  weight  minimiza- 


tion, speed  which  was  computed  in  step  3,  and 
number  of  parallel  paths  in  armature  winding,  the 
selection  of  terminal  voltage  will  decide  the 


number 
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To  determine  the  motor  performance  parameters,  t 
following  procedures  are  required: 

a.  Compute  the  diameter  of  armature  winding. 


b.  Compute  the  number  of  turns,  current,  resistance, 
and  inductance  of  armature  winding. 

c.  Compute  the  voltage  constant  and  torque  constant. 
The  computed  parameters  of  motor  inertia,  armature 
resistance,  armature  inductance,  voltage  constant, 
torque  constant,  and  damping  coefficient  will  be 
used  as  input  parameters  for  the  design  of  a control 
system  which  is  out  of  the  scope  of  this  disser- 


5.1  Formulation  of  Operation  Time  Enhancement 


5.1.1  Objective  Function 

We  consider  a single  motor  actuating  the  first  or  the 
second  joint  of  a multi-actuator  system.  We  presume  that  for 
the  goals  of  this  research,  all  parameters  of  the  system  may 
be  averaged.  Weight  minimization  of  the  actuator  at  each 
link  will  have  substantial  effect,  since  the  accelerating 


masses  will  be  reduced  correspondingly.  Thus,  minimizing 
the  weight  and  time  of  operation  for  a single  motor  should 
lead  to  a substantial  improvement  in  productivity  for  the 
whole  manipulator. 

Consider  the  motion  of  single  joint  of  actuator-ger- 
load  assembly  as  shown  in  Fig.  5-2.  There  are  two  basic 


t rotation 


Schematic  diagram  of  electric 


motions  in  this  single  joint — raising  arm  motion  and  waist 
motion.  To  raise  the  manipulator  arm  with  load,  the  motor 
torque  will  cope  with  the  weight  of  the  arm,  weight  of  the 
load,  and  acceleration  of  load  and  manipulator.  The  joint 
system  will  be  formulated  with  a rigid  body  structure,  and 
the  bending  capacity  of  gear  teeth  will  limit  the  applied 
motor  torque. 

According  to  the  D'Alembert's  principle,  the  mechanism 
must  be  in  equilibrium  under  the  influence  of  the  external 
and  internal  forces,  so  the  following  torque  equilibrium 
equations  can  be  obtained. 


T j - load  torque  at  load  side 
T = gravity  torque  of  manipulator 

= WtL  Cos  0 (in  waist  motion  T = 0)  (5.3) 


. = weight  of  manipulator 

I = load  torque  at  motor  side 
= (JjOjp  + Bj0m)/g l * Tw/gr 
. = load  inertia  including  manipulator  i 


B = damping  coefficient  at  motor  side 
Bj  = damping  coefficient  at 


0_  = angular  displacement  at  motor  side 
©2  - angular  displacement  at  load  side 

From  eq.  (5.2)  and  eq.  (5.4),  we  can  get  the  following 
expression : 

Tm  - Jeff°m  + Beff®m  + V»r  (5'51 


Je£g  = effective  inertia  at  motor  side 

" Jm  + V«£ 

= effective  damping  coefficient 

= Bm  + V*? 


If  we  assume  that  this  system  is  operating  from  a 
limited  source  of  power  and  it  is  desired  to  have  the  system 
change  from  one  state  to  another  in  minimum  time,  then  it  is 
necessary  at  all  times  to  utilize  all  the  power  available, 
that  is  to  use  '’bang-bang"  control  170]  . As  a result  of 
this  control,  the  magnitudes  of  acceleration  and  decelera- 
tion are  constant  and  maximum  values,  which  will  make  the 
velocity  trajectory  triangular.  From  the  triangular 
velocity  trajectory,  the  acceleration  torque,  damping  torque, 
load  torque,  and  motor  torque  can  be  determined  as  in  Fig. 
5-3.  As  can  be  seen  from  this  figure,  the  motor  torque 


f acceleration  torque,  damping  torque. 


gravity  torque  is  not  a simple  shape;  it  is  a function  of 


DA«"at  -A"’"  - X -t W ■ 


IV 
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>u,-Wvi|w 


Substituting  eg.  (5.7) 


.6)  yields 


(Qj/2)4  (Gj/2)6 

+ 216  9360  * ' 

- eg.  (5.9)  into  eg.  (5. 


The  acceleration  time  can  be  determined  by  solving  the 
guadratic  eguation  of  eg.  (5.10). 


The  operation 


3 two  times  the  acceleration  time  as 
symmetric  triangular  velocity  trajec- 


In  the  waist  motion,  the  gravity  torgue 
(5.6)  becomes 


-l  -I  * 


(5.13) 


the  integration  of  damping  torgue 
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=jfv 


■ Beffgc°l  (5.14) 

Substituting  eg.  (5.6)  and  eq.  (5.14)  into  eg.  (5.13)  yields 
the  following  operation  time. 


The  objective  function  for  raising  arm  motion  and  waist 
motion  are  formulated  as  eq.  (5.12)  and  eq.  (5. IS) . The 
constraints  in  this  time  reduction  problem  is  the  bending 
limit  of  gear  teeth.  In  this  study,  the  spur  gear  will  be 
used  for  the  computation  of  gear  bending  capacity. 

5.1.2  Bending  Capacity  of  Sour  Gear  Teeth  (711 

As  the  gear  rotates,  the  number  of  pairs  of  gear  teeth 
in  contact  simultaneously  varies  from  one  to  two  or  more. 
However,  the  torque  limit  can  be  determined  by  imposing  the 
most  unfavorable  conditions. 

Consider  the  ability  of  a single  tooth  to  carry  bending 
load  almost  at  the  end  of  tooth  as  shown  in  Pig.  5-4.  When 
the  force  P is  applied  at  point  A,  the  tangential  component 
of  F will  produce  a bending  moment  F.  1 at  the  base  of  the 

the  bending  limit  of  gear  tooth  will  be  determined  as  fol- 
lows. The  stress  of  a beam  can  be  expressed  as 
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when  the  gear  material,  diametral  pitch,  and  diameter 
of  pitch  circle  are  given,  the  bending  capacity  can  be  com- 


5.2  Time  Reduction 

Based  on  the  procedure  of  time  reduction  illustrated 
in  Fig.  5-1,  the  detailed  analysis  in  each  step  will  be 
discussed. 

5.2.1  Functional  Relationship 

The  functional  relationships  between  inertia  and  torque 
can  be  obtained  from  the  result  of  weight  optimization  in 
Chapter  IV.  The  relations  can  be  approximated  using  linear 
regression  of  the  SAS  package. 


Jm  = 0.00023Tr1,2364 

These  relations  express  the  dependencies  between  the  param- 
eters when  the  motor  is  in  optimum  configuration  with  respect 
to  minimum  weight  in  normal  and  "inside-out"  structure. 

Figure  5-5  shows  the  actual  data  and  approximated  data  of 
inertia  with  different  values  of  torque,  which  are  quite 
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: approximation  : actual  data 


Fig. 


Approximation  of  inertia  as  a function  of 
motor  ^ torque  and  actual  data  from  the  result 


5.2.2  Computation 


raising  arm  motion  is 


Minimize 


We  will  assume  that  the  effective  damping  coefficient  is 
zero)  then  eq.  (5.18)  of  raising  arm  motion  will  be  simpli- 


n bending  limit. 


Among  the  variables  in 
parametric  data  for  the 


[.  (5.19),  the  load  angular  displace- 
! given  as  input  parameters.  The 
•DMA  robot.  Unimate  500,  are  listed 


in  Table  5-1.  The  motor  inertia,  J , was  expressed  as  a 
function  of  motor  torque. 

When  the  constraint  is  not  applied  to  this  problem  to 
see  the  general  trend  of  time  function,  then  the  time  func- 
tion of  eq.  (5.19)  can  be  plotted  in  three  dimensions  with 
motor  torque  and  gear  ratio  as  variables  in  Fig.  5-6.  As 
can  be  seen  from  this  figure,  the  operation  time  of  the 
manipulator  is  reduced  as  the  motor  torque  becomes  larger. 


In  a given  motor  torque,  the  optimum  gear  ratio  can  be 
found  to  have  minimum  operation  time  from  this  figure.  The 
optimum  gear  ratio  can  be  determined  analytically  as 


3t0 
3g ” 


So  when  the  gear  ratio  in  each  motor  torque  is  deter- 
mined, then  the  operation  time  of  raising  arm  motion  can  be 
drawn  as  shown  in  Fig.  5-7.  We  can  see  from  this  figure 
that  the  operation  time  of  the  manipulator  using  the  "inside- 
out"  structure  motor  is  63%  less  than  the  manipulator  using 


normal  structure  motor.  When  the  maximum  bending  limit 
constraint  is  applied  to  this  figure,  the  maximum  available 
torque  is  the  optimum  torque  and  the  gear  ratio  will  be 
computed  successively. 

From  eq.  (5.17)  and  eq.  (5.18),  the  bending  limit 
torque  can  be  computed  with  different  values  of  diametral 
pitch  and  diameter  of  pitch  circle  as  shown  in  Fig.  5-8. 

When  the  diametral  pitch  of  the  gear  - 10  and  the  daimeter 
of  the  pitch  circle  D = 2.5  inches  (6.25  cm),  the  limit 


torque  is  determined  from  Fig. 
operation  time  is  determined  fr 
using  normal  structure  of  motor 


Pig. 


Operation 


-PW3 

= 7^ 
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3.  The  voltage  constant,  K , and  torque  constant,  K.  , 
will  be  computed  as 


5.3  Numerical  Example 

Using  the  Unimate  500,  the  design  of  a DC  PM  motor  in 
a single  joint  motion  is  computed  by  those  procedures  with 
a result  shown  in  Table  5-2.  The  "inside-out''  structure 
motor  has  better  characteristics  in  all  areas. 

Based  on  the  analyses  and  the  results  of  DC  PM  motor 
design  for  minimum  time  of  operation,  the  following  con- 
clusions can  be  obtained. 

1 . The  operation  time  is  reduced  by  matching  the  load 
and  the  actuator  conditions  with  maximum  bending 
capacity  of  gear  tooth,  and  the  problem  minimizing 

solved  via  hierarchical  design. 

2.  The  operation  time  of  raising  arm  motion  using 
"inside-out"  structure  of  motor  is  633  less  than 


design  is  developed  based 


the  requirements  of  TYPE  Ila  design. 


CHAPTER  VI 

MOTOR  DESIGN  STRUCTURES, 

CONCLUSIONS  AND  FUTURE  WORK 

Based  upon  the  results  of  analyses  in  this  dissertation, 
this  chapter  will  recommend  the  design  structures  of  a motor 
for  a multilink  manipulator.  The  conclusions  of  this  work 
are  presented  and  suggestions  for  future  work  are  discussed. 

6.1  Design  Structure 

Figure  6-1  describes  the  structure  of  motor  design  for 
the  actuators  in  a miltilink  manipulator.  Actually,  this 
design  structure  is  the  integration  of  discussions  in  Chap- 
ters IV  and  V.  In  other  words,  the  operation  time  is  reduced 
to  find  the  most  suitable  motor  torque  after  the  first  step 
optimization  of  motor  design  is  finished.  This  first  step 
design  determines  the  values  of  main  design  variables.  When 
the  time  reduction  is  finished,  then  the  second  step  of  motor 
design  is  done  to  find  the  detailed  motor  parameters  and 
performance  parameters  such  as  diameter,  current,  number 
of  turns  in  the  armature  winding,  resistance  and  inductance 

damping  coefficient  is  hard  to  compute,  so  this  may  be 
obtained  from  measurement  or  from  similar  machines. 
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6 . 2 Conclusions 


The  following  conclusions  may  be  obtained  from  this 
dissertation: 

1 . The  study  of  the  state  of  the  art  in  motor  design 
leads  to  the  hierarchical  optimum  design  structure 
which  divides  the  optimization  problem  into  several 
sub-optimization  problems.  This  decoupled  optimi- 
zation will  simplify  the  design  noticeably  and  has 
the  advantage  of  flexibility;  detailed  design  that 
may  be  different  from  different  purposes  can  be  done 
using  the  same  main  program. 

2.  The  notion  "optimum  configuration"  is  introduced 
which  is  defined  as  the  length/diameter  ratio  lead- 
ing to  the  minimization  of  the  weight  of  the  motor 
or  maximization  of  its  efficiency.  The  optimum 
configuration  is  computed  theoretically  for  a num- 
ber of  motors.  It  is  found  that  not  all  of  the 
industrially  manufactured  motors  satisfy  the 
criterion  of  optimum  configuration. 

3.  Prom  the  database  we  found  that  the  DC  PM  motor  has 
the  highest  T/W  ratio  among  any  other  type  of  motor 
as  shown  in  Table  6-1,  and  the  best  configuration 
of  DC  PM  motor  manufactured  in  industry  is  deter- 


ef f iciency 


optimum  configurations  which  may  be 
different  purposes,  and  the  optimum  configuration 
in  weight  minimization  is  not  changed  although  the 
power  rating  is  changed. 

The  result  of  optimum  design  of  the  DC  PM  motor 
shows  that  the  DC  PM  motor  has  much  better  character- 


istics than  the  wound  motor,  and  that  the  most 
desirable  motor  for  the  robotic  applications  or  for 


PM  "inside-out''  motor  with  SmCo  which  will  upgrade 
the  performance  characteristics  while  maintaining 
the  technical  requirements  of  the  actuator. 

The  operation  time  is  reduced  by  matching  load  and 
actuator  conditions  and  the  hierarchical  structure 
of  motor  design  is  especially  necessary  to  simplify 
the  problems  of  multiple-goal  optimization  and  to 
reduce  the  computing  time. 


1200  lines  of  Fortran  codes  using  the  nonlinear 
programming  algorithm,  pattern  search  and  Fibonacci 
search,  can  be  applied  to  any  other  types  of  machine 
by  just  modifying  the  usbprogram,  "DESIGN,"  which 


differs  for  different  types  of  motors. 

The  result  of  optimum  motor  design  can  be  applied 
to  the  actual  manufacturing  of  the  motor  and  will 
improve  the  motor  cost  effectiveness  and  other 


characteristics  while  increasing  the  productivity 
of  robotic  applications. 


rotor  dimensions 
we  can  suggest  a 

efficiency. 


for  Future  Research 

f the  database , the  dimensions  of 
outside  case  dimensions  which  are 
he  rotor  dimensions.  If  the  exact 
can  be  obtained  from  the  industry, 
more  accurate  optimum  configuration 
minimum  weight  and  maximum 


If  the  motor  is  actually  manufactured  from  the 
result  of  optimum  design  and  all  the  performance 
parameters  and  temperature  rise  are  measured  from 
this  manufactured  motor,  the  optimum  design  of  the 
motor  can  be  further  improved  by  comparing  the  com- 
puted parameter  with  the  measured  parameters. 

The  operation  time  was  reduced  by  matching  the  load 
and  actuator  conditions  with  the  concept  of  an 

motor  torque  will  provide  more  accurate  results. 


APPENDIX  A 
TRANSFORMER 


Figure  4-2  shows  the  simplified  geometry  of  C - core 
transformer.  If  the  flux  varies  sinusoidally  with  time  the 
primary  and  secondary  induced  voltages  are 


- rms  primary  and  secondary  voltages 

= number  of  primary  and  secondary  windings 
= frequency  of  power  source 

- peak  value  of  magnetic  flux 


The  window  area  which  is  the  winding  area  of  transformer 
is  assumed  to  be  divided  equally  between  primary  and 
secondary  windings,  therefore  the  total  number  of  turns  in 
the  primary  winding  can  be  determined  as 


k = window  utilization  factor 
D « diameter  of  wire 

In  order  to  calculate  the  current  in  any  winding,  we 
considered  the  current  density  at  first.  The  losses  of  a 
transformer  is  proportional  to  the  cube  of  geometrical 
dimension  and  the  square  of  current  density  [731. 
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The  area  product  which  is  an  important  variable 
can  be  expressed  as, 


for  relating 
of  windings, 


Ap=  (window  area)  x (cross 
- LcDm2(Dc-Dm)/2 


n be  expressed  a 

ss=  KV'75j2 


(A. 5) 


The  average  temperature  increase  of  a transformer  winding  is 
proportional  to  the  energy  losses  and  inversely  proportional 
to  the  surface  area  of  transformer.  As  a result  the 
temperature  rise  may  be  written  as. 


KlV'5 


(A. 6) 


Thus  at  a given  temperature  rise  the  current  density  is 
expressed  as  a function  of  area  product . 


induced  output 


:-ampere(VA) 


given  by 


and  the  losses  of  transformer  can  be  determined  as 
follows.  Mean  turn  length  of  each  winding  is 


n primary  and  secondary  windings  b 


(rlVV4) 

Resistance  of  each  winding  is  then 

R _ 4pLe (De+  Dm) <Dc~  Dm) ku  . 

OtDw2/4)2 

And  the  copper  losses  in  each  winding  can  be  expressed 
t _ pLc(Dc+  Dm)  (De-  1^)0 .75^^220 ,25 


m magnetic  flux  density  is  decided 
core  losses  per  unit  weight 


(A 


L ' . = core  losses  per  unit  weight 


When  the  permanent  magnet  is  used  as  a fj 
source  as  in  Pig.  4-5,  the  reluctance  of  the 
negligible  in  field  computation  for  this  kind  of  s 
Therefore  an  accurate  mathematical  model  in  field  determina- 
tion must  be  employed  to  take  all  the  parameters  into 

there  are  no  free  currents  inside  the  stator.  As  a result, 
the  flux  pattern  in  a steady  state  is  not  changed,  which 
enables  us  to  solve  the  Laplace  equation. 

The  governing  equations  for  this  case  are 


(B.l) 


Because  the  curl  o 
can  be  defined  as. 


, a scalar  magnetic  potential 


From  equations  (B.l)  and  (B.2)  with  B = uH,  the  magnetic 
scalar  potential  obeys  the  Laplace  equation 
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In  this  analysis  the  variation  of  scalar  potential  along  the 
z-axis  is  neglected.  And  if  only  the  fundamental  component 
of  harmonics  of  stator  winding  current  is  considered,  the 
three  phase  windings  approximates  a surface  current 
distribution  of  the  form 

K(8  t)  = Be  K a 5 (wt-m$)  ; . 


The  magnetic  potential  within  the  s 
sinusoidally  varying  with  time. 

f(r,t>  ■=  Re 

Substituting  eq.{B.6)  into  eq. (B.4) 
coordinates  yields 
X _d{,  dlirl,.  mf 


The  magnetic  scalar  potential  in  different  regions  such  as 
rotor,  airgap,  stator,  and  outside  the  motor  can  be 


expressed  accordingly. 
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Vr>  = V“  + Vr"  ! ais rsa2 

Vr)  = A4rra  + A5/rm  : 


?4(r)  = Ag/rm  ' I!a3  (B  . 11) 


«!♦=  H2  ' Blr”  B2r 


H2*'  Hl*=  V B2r“  B3r  (B-13> 

H3*-  H4*  ' B3r=  B4r  <B‘14> 


(jXlV“  + e3  (B . 15) 

= [jKz(X2rm+^)+e:iS(XJrm+^!]e:ilwt"'n41  (B.16) 

= DKZ *X4 rI”+^ml  +e35  (X4rm+^|)]ej  (B.17) 

= !jKs,X6/r,n  + e35X£/rn,)e3(wt'rai|’>  IB. 18) 


v *sv^ 

>s" 


V k-x3 


X2  * *J*i  + *3 


l'3/a*m 


■ 5 [a*"x2  (1-Wo/v)  +(1+V„/M)  Xj 

• 5 [■a"xJ ( +<  l*Wm  X 


UoM/  <M-Uo)ma“' 


-(U+U0!/  (Wl>o>max 


ax=  d/2 


Among  those  magnetic  fields  only  the  radial  component  of 
magnetic  field  in  the  airgap,  HJr,  and  the  tangential 
component  of  magnetic  field  in  the  stator  are  of  interest  to 
us.  The  magnetic  flux  in  the  stator  is  obtained  from  the 

*3ru.fXH3tdr  dL 

=ML{X5  < 1/a 1/a”)  -XJ  (a“-a“)  I^CosIwt-i#) 

+ X^ ( 1/ a^-l/ a™) “X^ (a™- a”)  Sin(wt-m$-6) ) (B.23J 

The  magnetic  flux  linkage  that  is  cut  by  the  stator  winding 


The  coenergy  stored  in  the  airgap  is 


The  torque  exerted  on  the  PM  rotor  becomes 

The  pullout  torque  and  the  average  torque  of  a three  phase 
m pole  pair  motor  are  then 

Tp0=  3 . 332mvi„ LKw  Hi  (B.27) 

Tav=  -|i0LKwNi(XJ-XJrra)  Sin  6 (B.28) 

When  the  number  of  pole  is  selected  the  angular  speed  of 
rotor  in  a synchronous  motor  is  determined  and  the  average 
output  power  can  be  computed  accordingly  as. 


APPENDIX  C 

PLOW  CHART  OF  COMPUTER  PROGRAM 
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''  PATTRN’ 
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STOMJflBS  (OBJX-OBJ)  /C8J 


C=D 


"THERML" 


F COMPUTER  PROGRAMS 


For  the  execution  of  this  computer  program  the  Harris 
computer  program  are  as  follows: 

a.  Turn  on  the  terminal. 

b.  Hit  <control-G>  Prompt  (USER  i ?) 

c.  Enter  user  name  <spacexcontrol-Z>  password  <return> 
2.  Compiling  the  computer  program  "OPTIM” 


b.  SAUF77  OPTIM  <return> 

c.  VU.R  <return> 


Begin  <return> 

Assigning  each  file  to  ea 
1S=0PTDAT  <return> 


d.  AS  4=EVE  <return> 

e.  AS  6=*  <return> 

f.  AS  7=STEP  <return> 


Execution 
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